Convex Hull of Planar H-Polyhedra 



Axel Simon and Andy King 
Computing Laboratory 
University of Kent, CT2 7NF, UK 
{a.simon, a.m. king}@ukc. ac.uk 

February 1, 2008 



Abstract 

Suppose (Ai,Si) are planar (convex) H-polyhedra, that is, Ai G R' l * x2 
and Ci £ R n * . Let Pi — {x £ R 2 \ AiX < Ci} and n = n\ + n^. We present 
an 0(n log n) algorithm for calculating an H-polyhedron {A, c) with the 
smallest P = {x G R 2 | Ax < c\ such that Pi U P 2 C P. 

Keywords: convex hull, computational geometry 

C. R. Categories: 1.3.5 [Computational Geometry and Object Modeling]: 
Boundary representations; Geometric algorithms, languages, and systems, 1.3.6 
[Methodology and Techniques]: Graphics data structures and data types, F.3.1 
[Specifying and Verifying and Reasoning about Programs]: Invariants, Mechan- 
ical verification. 



1 Introduction 

The convex hull problem is classically stated as the problem of computing the 
minimal convex region that contains n distinct points {(xi, ?/i)}™ =1 in the Eu- 
clidean plane K 2 . The seminal work of Graham showed that the convex 
hull problem can be solved in O(nlogn) worse-case running time. It inspired 
many to elaborate on, for example, the three and more dimensional case, spe- 
cialised algorithms for polygons, on-line variants, etc. [51 HUj. The convex hull 
of polytopes (bounded polyhedra) can be calculated straightforwardly by taking 
the convex hull of their extreme points. However, calculating the convex hull 
for polyhedra turns out to be more subtle due to a large number of geometric 
configurations. Even for planar polyhedra, the introduction of rays makes it 
necessary to handle polyhedra such as a single half-space, a single ray, a single 
line, two facing (not coinciding) half-spaces, etc., all of which require special 
handling in a point-based algorithm. The problem is exacerbated by the num- 
ber of ways these special polyhedra can be combined. In contrast, we present 
a direct reduction of the convex hull problem of planar polyhedra to the con- 
vex hull problem for a set of points 0]. By confining all input points to a box 



and applying the rays to translate these points outside the box, a linear pass 
around the convex hull of all these points is sufficient to determine the result- 
ing polyhedron. By adopting the classic Graham scan algorithm, our algorithm 
inherits its 0(n log n) time complexity. The standard tactic for calculating the 
convex hull of TJ-polyhedra is to convert the input into an intermediate ray and 
vertex representation. Two common approaches to this conversion problem are 
the double description method [7] (also known as the Chernikova algorithm [3]) 
and the vertex enumeration algorithm of Avis and Fukuda j^j ■ The Chernikova 
method leads to a cubic time solution for calculating the convex hull of planar 
iJ-polyhedra |S] whereas the Avis and Fukuda approach runs in quadratic time. 

The remaining sections are organised as follows: A self-contained overview 
of the algorithm, together with a worked example, is given in the next Section. 
A formal proof of correctness is given in Section [3J Section 0] concludes. 

2 Planar Convex Hull Algorithm 

The planar convex hull algorithm takes as input two -ff-polyhedra and outputs 
the smallest if-polyhedron which includes the input. The ^/-representation of 
a planar polyhedron corresponds to a set of inequalities each of which takes 
the form ax + by < c, where a,b,c € R and either a ^ or b ^ 0. Let Lin 
denote the set of all such inequalities. The vector (a, b) is orthogonal to the 
boundary of the halfspace induced by ax + by < c and points away from the 
feasible space. This vector induces an ordering on halfspaces via the orientation 
mapping 9. This map 9 : Lin — > [0, 2n) is defined such that 9(ax + by < 
c) = tfi where cos(V') = aj\Ja 2 + b 2 and sin(^) = bj\/a 2 + b 2 . The mapping 9 
corresponds to the counter-clockwise angle which the half-space of x < has 
to be turned through to coincide with that of ax + by < c. Sorting half-spaces 
by angle is the key to efficiency in our algorithm. However, 9 is only used for 
comparing the orientation of two half-spaces. To aid the explanation of the 
algorithm, the concept of angular difference eiZe2 between two inequalities e\ 
and e-i is introduced as the counter-clockwise angle between 9{e\) and 9{e2)- 
More precisely eiZe2 = (9(e2) — 9(ei)) mod 2ir. Note that this comparator can 
be realized without recourse to trigonometric functions 9 . 

The algorithm makes use of a number of simple auxiliary functions. The 
function intersect{a\x + b\y < c\, a2X + biy < c-i) calculates the set of intersec- 
tion points of the two lines a\x + b\y — c\ and a2X + b^y = c%. In practice an 
implementation of this function only needs to be partial since it is only applied 
in the algorithm when the result set contains a single point. The remaining 
auxiliaries are listed in Figure ^ The connect function generates an inequal- 
ity from two points subject to the following constraints: the halfspace induced 
by connect(pi,p2) has pi and pi on its boundary and if p\,p2,py, are ordered 
counter-clockwise then p$ is in the feasible space. The notation pT7P~2 is used 
to abbreviate connect{p\,p2)- Furthermore, the predicate saturates (p,e) holds 
whenever the point p is on the boundary of the halfspace defined by the inequal- 
ity e. Finally, the predicate inBox(s 7 p) determines whether the point p occurs 



1 function extreme(E) begin 
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ifE = {ax + by < c} then R := {(-a/Va 2 + 6 2 , -b/y/a 2 + b 2 }}; 


■5 


for ie [0, n — 1] do let ei = ax + by < c in begin 
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// are the intersection points of this inequality degenerated? 


7 


d pre := {0{ei) - 0(e»_i mo d n)) mod 2ir > tt V n = 1; 


8 


rfpost := (0(e,+i mod n) - 0(e*)) mod 2tt > ir V n = 1; 
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if d pre thcni? := i?U{(&/Va 2 + 6 2 ,-a/Va 2 + & 2 )}; 


10 


if d post then J? := i? U {(-o/Va 2 + b 2 , a/Va 2 + b 2 }}; 


11 


else V :=VU intersect^, e (i+ i) mod „); 


12 


if d pre A d pos t then := U {«} where u € {(x, y) | ax + 6y = 


13 


end 


14 


return (V, i?) 


15 end 


16 





17 function connect({x±, y±), (x2,y2)) 

18 return (y 2 - yi)x + (xi - x 2 )y < (y 2 - Vi)xi + (%i - x 2 )yi 
19 



20 function saturates{{x\, yi), ax + by < c) 

21 return (axi + byi = c) 
22 

23 function inBox(s, (x,y)) 

24 return |x| < s A \y\ < s 



Figure 1: Convex hull algorithm for planar polyhedra 



within a square of width 2s that is centred on the origin. 

The algorithm divides into a decomposition and a reconstruction phase. 
The hull function decomposes the input polyhedra into their corresponding 
ray and vertex representations by calling the function extreme in lines 3 and 
4. The remainder of the hull function reconstructs a set of inequalities whose 
halfspaces enclose both sets of rays and points. The functions extreme and 
hull are presented in Figures ^ and [21 respectively. The algorithm requires the 
input polyhedra to be non-redundant. This means that no proper subset of 
the inequalities induces the same space as the original set of inequalities. The 
algorithm itself produces a non-redundant system. 

To illustrate the algorithm consider Figured! The polyhedron E = {eo, ei, e 2 } 
and the polytope E' = {e' , . . . ,e' 5 } constitute the input to the hull function. 
They are passed to the function extreme at line 28 and 29. Within extreme the 
inequalities of each polyhedron are sorted at line 2. Note that for ease of pre- 
sentation the indices coincide with the angular ordering. The loop at lines 5-13 
examines the relationship of each inequality with its two angular neighbours. If 
d pos t is false, the intersection point intersect^, e^ + i) mo d n ) is a vertex which is 



25 function hull(Ei, E 2 ) begin 




26 


/ / assertion: each Pj is satisfiable and non-redundant 




27 


if Ei = V P 2 = then return 0; 




28 


(Pi,Ri) := extreme(Ei); 




29 


(P 2 ,R 2 ) ■= extreme(E 2 ); 




30 


P := PiUP 2 ; 




31 


R := RiUR 2 ; // Note: \R\ < 8 




32 


s := max{\x\,\y\ (x,y) G P} + 1; 




33 


/ / add a point along the ray, that goes through x, y 




34 


II and is outside the box 




35 


Q ■= P; 




36 


for (x,y,a,b) e P x R do Q := QU{(x + 2s/2sa, y + 2\/2s6)}; 




37 


/ / construct four inequalities in the zero dimensional case 




38 


if Q = {(x\, j/i )} then return {x < x\,y < yi, —x < —Xi, —y < 


-yih 


39 


/ / the centre of gravity q p is feasible but not a vertex (since \Q\ > 1) 


40 


% : = CE( x ,y)eQ x /\Q\^( x , y )eQy/\Q\^ 




41 


II q p is pivot point for sorting: Mi G [0, n— 2] . 9(q p ,qi) < 0(q p , 


Qi+i) 


42 


(go, • • -,qn-i) ■= sort(g p ,Q) 




43 


/ / identify the m vertices q^ where < k < . . . < fc m _i < n 




44 


(Qk ,- ■ ■ ,9fe ro _i) := scan((g , • • ■ ,q n -i}) 




45 


Pros := 0; 




46 


for i € [0, m — 1] do begin 




47 


let (xi,t/i) =g fei , (x 2 ,t/ 2 ) =qk (i+1)modm 




48 


let e= connect ((x i, y±), (x 2 , y 2 )) 




49 


If add e to P ros if q kz or g fe( . +1) mod m is in the box... 




50 


add := inBox{s, (x\,y\)) V inBox(s, (x 2 ,y 2 )) V to = 2; 




51 


j := (fcj + 1) mod n; 




52 


while ^arfrf A j 7^ fc i+ i do begin 




53 


// ...or any boundary point is in the box 




54 


add := saturates(qj, e) A inBox(s, qj); 




55 


j ■= (j + 1) raod n; 




56 


end; 




57 


if to = 2 A inBox(s, (x\, y\)) then 




58 


if yi = 2/2 then E Tes :=E ICS U {s5n(xi — x 2 )x < sgn(xi — 


x 2 )x x } 


59 


else P ros :=P ros U {s5?i(yi - j/ 2 )y < sgn(y 1 - 


2/2)2/1} 


60 


if ac?c? then P ros := P ros U {e}; 




61 


end 




62 


end 




63 


return P res 




64 end 





Figure 2: Convex hull algorithm for planar polyhcdra 
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polyhedron. 



Figure 3: The different stages of the polyhedra convex hull algorithm. 



added at line 11. Conversely, if d post is true, the intersection point is degenerate, 
that is, either E contains a single inequality or the angular difference between 
the current inequality and its successor is greater or equal to 7T. In the example 
two vertices are created for E, namely v\ and V2 where {vi} = intersect(eo,ei) 
and {112} = intersect(ei,e2)- The intersection point intersect^, eo) is degener- 
ate, thus it is not added to V; in fact the point lies outside the feasible space. Six 
vertices are created for E 1 . Rays are created at line 9 and 10 if the intersection 
point is degenerate. The two rays along the boundaries of e^ and eu+i) mod n 
are generated in loop iteration i when d pos t is true and iteration (i + 1) mod n 
when d pre is true. In our example d pos t is true for e%, generating a ray along 
the boundary of e2 which recedes in the direction of the first quadrant, whereas 
d pre is only true for eo yielding a ray along eo which recedes towards the second 
quadrant. No rays are created for the polytope E' . 

In general both flags might be true. In this circumstance the current inequal- 
ity ei cannot define a vertex. In this case an arbitrary point on the boundary 
of the halfspace of a is created at line 12 to fix its representing rays in space. 
Another case not encountered in this example arises when the polyhedron con- 
sists of a single halfspace (\E\ = 1). In this case a third ray is created (line 4) to 
indicate on which side the feasible space lies. Note that the maximum number 
of elements in R never exceeds four, which occurs when the input defines two 
facing halfspaces. 

The remainder of the hull function is dedicated to the reconstruction phase. 
The point and ray sets, returned by extreme, are merged at line 30 and 31. At 
line 32 the size of a square is calculated which includes all points in P. The 
square has (s, s), (— s,s), (s,— s), (— s, — s) as its corners. The square in the 
running example is depicted in all three frames of Figure 01 and the origin is 
marked with a cross. Each point p G P is then translated by each ray r G R 
yielding the point set Q. Translated points appear outside the square since all 
normalised rays are translated by the length of the diagonal 2y2s of the square. 
The translation process for the worked example is depicted in the second frame. 
Line 38 is not relevant to this example as it traps the case when the output 
polyhedron consists of a single point. Line 40 calculates a feasible point q p of 
the convex hull of Q which is not a vertex. This point serves as the pivot point 
in the classic Graham scan. First, the point set Q is sorted counter-clockwise 
with respect to q p . Second, interior points are removed, yielding the indices of 
all vertices, in the case of the example ko, . . . , k-j. What follows is a round-trip 
around the hull which translates pairs of adjacent vertices into inequalities by 
calling connect at line 48. Whether this inequality actually appears in the result 
depends on the state of the add flag. In our particular example the add flag is 
only set at line 50. Whenever it is set, it is because one of the two vertices lies 
within the square. The resulting polyhedron consists of the inequalities q^ 2 , , 
qk 3 , q^ i and qjeliOks which is a correct solution for this example. 

In general, the reconstruction phase has to consider certain anomalies that 
mainly arise in outputs of lower dimensionality. One subtlety in the two dimen- 
sional case is the handling of polyhedra which contain lines. This is illustrated 
in Figure0]where the two inequalities eo, e\ are equivalent to one equation which 



Figure 4: The resulting convex hull contains a line. 



defines a space that is a line. Observe from the second frame that no point in the 
square is a vertex in the hull of Q. Therefore the predicate inBox does not hold 
for the two vertices qk 2 and qk 3 and the desired inequality qkT7<lkI is n °t emitted. 
Similarly for the vertices and qk - However, in such cases there always ex- 
ists a point in p 6 Q with q p , q^/Lqp^p < q p , qj^£%, Qk (i+1) mod ~ which lies in the 
square. Thus it is sufficient to search for an index j € [fej + 1, fc(,+i) mo d m ~ 1] 
such that qj is both in the square and on the line connecting the vertices qt t 
and qk (i+1) modm - The inner loop at lines 52-56 tests if Q contains such a point 
and sets add appropriately. 

The one dimensional case is handled by the m = 2 tests at line 50 and 57. 
Figure [U illustrates why the test in line 50 is necessary. Suppose E\ and are 
given such that extreme(Ei) = ({94, 95}, 0) and extremeiE^) = ({93}) { r i ~ r }) 
where r is any ray parallel to the line. Observe that all points are collinear, thus 
the pivot point is on the line and a stable sort could return the ordering depicted 
in the figure. The correct inequalities for this example are E reB = {qo, qs, q~g, qo}. 
The Graham scan will identify qk = qo and qk t = q& as vertices. Since there 
exists j € [ko + 1, fci — 1] such that inBox(s, qj) holds, qo, q$ G -E TO s- In contrast, 
although there are boundary points between qs and qo the loop cannot locate 



Figure 5: Handling the one-dimensional case. 



them due to the ordering. In this case the m = 2 test ensures that add is set, 
guaranteeing that q$, qo € E rcs . 

The output polyhedron must include q ki as a vertex whenever inBox(s, q^) 
holds. If inBox{s, q ki ) holds, the algorithm generates ej_i = qk (i _ 1} mod ~ , <7fci and 
e » = Qki,qk (i+1) modm - If e-i-i^-Zi < it, then {q ki } = intersectfa-i, ej) and the 
vertex is realized. Observe that if m = 2, ei_iZe.; = 7r which necessitates an 
additional inequality to define the vertex q ki . This is the role of the inequality 
generated on lines 58 or 59. Observe that this inequality e guarantees ei_iZe < 
7r and e/-&i < tt which is sufficient to define q ki . 

The zero dimensional case corresponds to the case of when both input poly- 
hedra consist of the single point v. Line 38 traps this case and returns a set 
of inequalities describing {v}. Observe that the zero and one dimensional case 
only require minute changes to the general two dimensional case. 

As a note on implementation, observe that the search for a pivot point at 
line 40 can be refined. One method for finding a definite vertex is to search for 
a point with extremal coordinates [J. However, this process requires all points 
to be examined. The presented algorithm follows Graham ;4 in creating an 
interior point as the pivot point. This does not necessarily require the whole 
point set to be examined. By choosing two arbitrary points qi,q2, it is sufficient 
to search the point set for a point qi which does not saturate the line qT7Q2- The 
centre of the triangle q\,q%, q% is guaranteed to be an interior point of Q. 

Note also that the sorting in extreme is unnecessary if the input inequalities 
are consecutive in terms of angle. In particular, the output of one run of the 
algorithm can serve as an input to another without applying the sort at line 2. 

Finally observe that the inner loop at lines 52-56 can often be skipped: if the 
line between q ki and qk {i+1} mod m does not intersect with the square, inBox(s, q) 
cannot hold for any q £ Q. Hence add will not be set at line 54 and the inner 
loop has no effect. 



3 Proof of Correctness 



Section 13. II introduces the mathematical language necessary for expressing the 
two parts of the proof. The proof itself reflects the structure of the algorithm: 
Section 13.21 concerns the conversion of planar ii-polyhedra into their ray and 
point representations; Section argues that the reconstructed polyhedron en- 
closes the points and rays generated from the two input polyhedra, and yet is 
also minimal. 

3.1 Preliminaries 

A polyhedron is a set P C R d such that P = {x £ R d | Ax < c} for some 
matrix A 6 R nxd and vector c £ R n . An H-polyhedron is a pair (A, c) where 
A e R nxd and c £ E™ that is interpreted by [.] as the polyhedron {(A,c)j = 
{x £ R d I Ax < <?}. For brevity we manipulate iZ-polyhedra as a finite set of 
inequalities E = {a\ ■ x < Ci, . . . , a n ■ x < c n } which is equivalent to the matrix 
representation (A, c) with A = (oj., . . . , a„) T and c = (c\, . . . , c n ) T . E is said to 
be satisfiable if [E] ^ and non-redundant if Ve £ E . \E \ {e}] ^ \E\. Two 
inequalities e, e! coincide, written c(e, e'), iff there exists a € R d ,c € M. with 
[{e,e'}l = {x G R d | a- x = c}. 

The convex hull of a finite set of points P = {pi, . . . ,p n } CM. d is defined as 
conv(P) = {J27=i ^iPi I — ^i^Y^i=i = !}■ Moreover, the cone of a finite set 
of vectors R — {fi, . . . , f m } C R d is defined as cone(R) — {YliLi ^i^i I < Aj}. 
The Minkowski sum of two sets X,Y C M. d is defined a,sX + Y = {x + y\ 
x e X A y e Y}. Let vert(S) = {p £ S \ p <£ conv(S \ {p})}, ray(S) = {r £ 
R d \{6} | S+cone{{r}) C S} and line(S) = {r £ ray(S) \ -r £ ray(S)} denote 
the vertices, rays and lines of a convex set S. The following result, accredited 
to Motzkin ^I] , relates the two classic representation of polyhedra. 

Theorem 3.1 The following statements are equivalent for any SCR 11 : 
1.5= conv(P) + cone(R) for some finite P, R C R d ; 
2. S = {x £ R d | Ax < c} for some matrix A £ R nxd and vector c£ R n . 

Our algorithm converts the two (planar) input i/-polyhedra Ei into their 
rays Ri and points Pi and calculates an 7J-polyhedron E with the smallest [E] 
such that [Ej D (E^ U \E 2 ]. In fact \E] = conv(P 1 U P 2 ) + cone(Ri U R 2 ). 

3.2 Decomposition 

The discussion of the function extreme is organised by the dimension of the 
polyhedron. Reoccurring or self-contained arguments are factored out in the 
following lemmata. The first lemma states how redundancy can follow from 
the angular relationship between three inequalities. Since extreme requires its 
input to be non-redundant, useful angular properties of the input inequalities 
flow from the lemma. 



Lemma 3.1 Suppose ej = aix + biy < a, i — 1, 2, 3. Let c\ = C3 = 0, C2 > 0. 
Tften [{ei,e 3 }] C [{e 2 }] if = 9{a) < 0(e 2 ) < 0(e 3 ) < n. 

Proof. To show |{ei,e3}] C |{e2}] it is sufficient to show [{ei,e3}] C \{a 2 x + 
b 2 y < 0}], thus let c 2 = 0. Furthermore, w.l.o.g. let e± = x < (since 6(e\) = 0), 
e 2 = a 2 x + y < (since 9(e 2 ) < tt) and e3 = 03^ + y < (since 0(e3) < 7r). 
Note that 02 = Aiai + A3<i3 and b 2 — \±bi + A3&3 has the solution Ai = a 2 — 03 
and A3 = 1. Due to b 2 = 1, a 2 = cot^ 1 (8(e 2 )) and similarly 03 — cot _1 (#(e3)). 
Because 0(e 2 ) < Q(e$) and cot is an anti-monotone on (0, 7r), it follows that 
a 2 > 03, hence Ai > 0. Let (x, y) satisfy e\ = a; < and e3 = 03a; + y < 0. Due 
to Aix < 0, Ai(x) + A3(d3iz; + y) = (Ai + 03)0; + y < = e 2 , thus e 2 holds, hence 
[{ei,e 3 }] C [{e 2 }] as required. ■ 

The following lemma states that there is an injection between vertices and in- 
equalities. This result is used to show that the loop in extreme does indeed 
generate all vertices of the input polyhedron. 

Lemma 3.2 Let E = {eo, . . . , e n —i} be non-redundant and ordered by 9 and 
v £ vertd^]). Then there exists e$ G E such that {v} — intersect^, e^) and 
e^X-ey < 7r where i' = (i + 1) mod n. 

Proof. Let E 1 C E contain those inequalities that v saturates. Note that 
\E'\ > 2, in particular there exist e,e' G E' with eZe' ^ {0,7r}, otherwise 
IE'] \ {v} is not convex. Choose e^e^ G E' such that e^e/,. < tt. Suppose for 
the sake of a contradiction that k ^ (i + 1) mod n. Then e 3 £ E \ E' exists 
with eiZej < ej/e&. W.l.o.g. assume that 0(ei) = and v = (0,0). Then 
ei^ej < eiZe k reduces to = 0(ej) < 9(ej) < 9(e k ) < tt. Lemma l3~Tl implies 
IE \ {ej}} C which contradicts the assumption that E has no redundant 
inequalities. Thus k = (i + 1) mod n = %' I 

The following result also builds on Lemma 13.11 and complements the previous 
lemma in that it states when adjacent inequalities have feasible intersections 
points. The lemma is used to show that extreme only generates points in \E\. 

Lemma 3.3 Let E = {eo, . . . , e„_i} be satisfiable, non-redundant and ordered 
by 9. For any e t £ E if ei^e (i+1 ) modn < tt then intersect (a, e (l+1) mod „) C 

Proof. Let e.; £ E. Since ei is not redundant in E the boundary of Ci intersects 
with the non-empty convex body \E \ {ei\\. Let e m £ E \ {a} such that 
^ S = intersected, e m ) C \E \ {e,}]. Note that if \S\ > 1 then E = {e i; e m } 
with eiZe m = tt, hence ejZe^+i) mo( j „ < 7r never holds. Let {v} = S. It remains 
to show that e, and e m are adjacent. Suppose that Ci/Le m < tt (e m Zei < tt is 
analogous). Assume for the sake of a contradiction there exists e; G E \ {e{\ 
such that ei^ei < eiZe m . W.l.o.g. v — (0,0) and 9(ei) = 0, hence = #(ej) < 
9(ei) < 9(e m ) < tt. Since v £ \E \ {ej}J, w G [e;] and thus c; > where 



e; = aix + biy < q. By Lemma l3.1l e; is redundant in E which is a contradiction. 
It follows that m— (i + 1) mod n. ■ 

While the previous lemmata concern points, the following lemma is a statement 
about rays. In particular, it states when inequalities give rise to rays. 

Lemma 3.4 Let E — {eo = a$x + boy < cq, . . . ,e„_i} be satisfiable, non- 
redundant and ordered by 9. Let i' = (i + 1) mod n, S\ = cone({(— bi, ai), 
(by,— ay)}) and S2 = ray(\E\). If eiLey > it or \E\ > 3 A e,Ze,v = ir then 
Si = 82- 

Proof. To show 5i C 82- Choose p € {EJ such that p saturates ej. For the 
sake of a contradiction assume = (—bi,cii) ^ ray ( [i?] ) . Thus there exists 
A > with p + Xri £ [Ej. Set A max = max{A > | p + Xr, E [E]}. Note 
that v = p + A max ri is a vertex and by Lemma 13.21 there exists e k € E with 
{v} = intersect(e k ,e( k+1 ) modn ) and e fe /e( fc+ i) mod „ < n. The latter implies 
that k =/= i, therefore i = (fc + 1) mod n. Hence e(,-_;n mo d n = does not 
contain the ray. W.l.o.g. let 0(e,) = and v = (0,0). Then = (0,1) and 
7T < 0(e(»_i) mod „) < 2tt, hence e (i _ 1} mod „ = a'x + b'y < with b' < 0. Thus 
a'x + b'(X + y) < 0, hence p + Ar; e |{e(j_i) mo d „}] which is a contradiction. 
Analogously for ry = (by,— ay). Now to show £2 C S^. Assume there exists 
re ray([Ej)\ cone({{-bi,ai),{bi',-ai>)}). Consider e l Le v > n. Since 0(ej) ^ 
9(ei') there exists Ai, A2 G K with r = Aj (— 6j, ai)+X2(bi' ,—ay). Assume A2 < 0. 
W.l.o.g. let {p} = intersect(ei,ei') = {(0,0)} and 9{ei) = thus let ei = x < 0. 
It follows that 7T < 9(ei>) < 2ir, by < and we set = a'^x — y < 0. Thus 
r = Ai<0, 1) + A 2 (-l, -4) = (-A 2 , Ai - X 2 a! i> ). Since A 2 < 0, p + r $ [{ej}] 
which contradicts r € rai/([_E]). Analogously for Ai < with 9(ei>) — 0. Since 
Ai > and A2 > 0, r G cone({(— 6i, ai), (bii , — ay)}) which is a contradiction. 
Now suppose eiZey — ir and \E\ > 3. W.l.o.g. let 9(e{) — 0, a = x < \u\ and 
ey = —x < 0. Thus Si — cone({(0, 1)}). Let (x r ,y r ) = r and p € [{e^, ej'}]. 
Observe that x r = otherwise there exists A > with p+ Xr ^ [{ej, ej'}]. Now 
assume y r < 0. Since \E\ > 3, there exists ej £ E such that tt < 9(ej) < 2ir and 
thus ej = ajX + bjy < Cj with bj < 0. Let p € [{ej}], then there exists A > 
such that p + Xr ^ [{e^}] which is a contradiction. ■ 

The correctness of the first stage of our algorithm is summarised by the 
following proposition. Note that the Minkowski sum conv(V) + always defines 
the empty space rather than the polyhedron conv(V). Thus a null ray (0,0) 
is required to represent a bounded polyhedron. The function extreme avoids 
adding this null ray for the sake of improved efficiency. In hull the translation 
of points at line 36 by the null ray is replaced by a simple copying step at line 
35. However, the null ray still manifests itself in the correctness results. 

Proposition 3.1 Let E C Lin be non-empty, finite, satisfiable and non-redun- 
dant. Then extreme(E) = (P, R) and {Ej = conv(P) + cone(R U {(0, 0)}). 



Note that the following proof handles polyhedra that contain lines as a spe- 
cial case. This distinction is not artificial, in fact Klee |S] observed that a closed 
convex set that does not contain lines is generated by its vertices and its extreme 
rays. (Extreme rays are rays that cannot be expressed by a linear combination 
of others.) In order to describe polyhedra which contain lines it is necessary to 
create points that are not vertices and rays which are not extreme. 
Proof. Let {eo, . . . ,e n -i} = E such that 0(eo) < 0{e\) < • • ■ < Q(e n -i)- Such 
an ordering exists since if 9(ei) = 9(ej) for some i ^ j then either [aj C [ej], 
thus [E\ {ej}} = [EJ or vice-versa which contradicts that E is non-redundant. 
Let (P,R) = extreme(E). 

• Suppose line(\E\) ^ 0. Note that if there exist ej, ej S E with e^/e^ < 7r 
then ime([£]) = 0, hence E = {e} or E = {e , ei} A e Zei = n. 

— Consider E = {e}. The loop is executed only once with d pre = 
dpost — true, thus the boundary point p and two rays rx,r% are 
added in line 12, 9 and 10, respectively. The boundary of \EJ is 
{p + AiTi + A2T2 | Xi > 0}. The ray r% added in line 4 points into the 
halfspace, thus [EJ = {p} + cone({n, r 2 , f* 3 }). 

— Consider E — {eo, e±} A eo/ei = n. For each e^ d pre — d pos t — true, 
hence for each e^ the loop generates two rays r, , r[ along the boundary 
of ej (lines 9, 10) and one boundary point Pi (line 11). The set 
conv({pi,p2]) is included in [EJ. Note that the rays generated in 
the second iteration are collinear to those in the first. Thus [£] = 
conv({pi,p 2 }) + cane({ri,r'i}). 

• Suppose E is zero dimensional, hence {Ej — {v}. Let Qfj = ejZe(j_)_i) mo d n . 

Since Y^i=o a i — 2tt, it follows that n > 3 since otherwise there exists 

1 such that ai > it and by Lemma I3.4I it follows that ray([_E]) ^ 0. 
By Lemma I3.ll it follows that 014 + ct(i+\) modn ^ which for all ej is 

Yh=0 ( a i + a (i+l) mod n) > n>TT. However, 717T < E"=0 ( a i + a (i+l) mod n) = 

2 Ya=o ai = ^ 7r ' hence n < 4. 

— Consider E — {eo, ei, ea, 63}. Observe that for all < i < 3, 
a, + a(j+i) mod n = 7T and intersect^, e( i+1 ) mod n ) = {v}, hence 
c ( e o,e2) A c(ei,e3). In all loop iterations d pr - e = d pos t = false thus 
only vertex {v} = P is generated (line 11). 

— Consider J5 = {eo, e±, ei\. Observe that for all i, eiZe( i+1 ) mo & n < ir, 
otherwise by Lemma 13.41 ray([i?]) ^ 0. Thus in each iteration i, 
dpost = d pre = false and P = {v} = intersect^, e^ l+1 ) modn ) is 
generated. Hence \E\ = {v} + cone({(0, 0)}). 

• Suppose E is one dimensional and line{\E\) = 0. Since the boundary of 
[EJ contains infinitely many points and \E\ = n (which is finite), there 
exists e € E such that p\,pi £ [EJ where p\ ^ pi and p\ and P2 saturate 
e. In fact there are infinitely many boundary points on the line between 



Pi andp2 which saturate e. Observe [EJ contains no interior points, hence 
[EJ consists only of boundary points. Therefore there exists e' £ E that 
saturates infinitely many of these points and for which c(e, e') holds. As 
we assume that line(\E\) = 0, \E\ < 2. Due to non-redundancy, there 
exists at most one ej € E with < eZe^ < it. Similarly for e'. Hence 
3 < \E\ < 4 follows. 

— Consider E — {eo,ei,e2}. Let i S [0,2] such that c(ei,e(i + i) modn) 
holds. On iteration i the loop generates a ray r along the boundary of 
[ej] in line 10. A collinear ray is generated in iteration (i + 1) mod n 
for e( i+1 ) mod n on line 9. It is in this and iteration (i+2) mod n where 
{v} = intersect(e {l+l) m0(i „,e (l+2 ) mo dn) = intersect(e {l+2) modn ,ei) 
is added to P. Thus [E] = {v} + cone({r}). 

— Consider = {eo, ei, e2, e3}. Let i € [0, 3] such that c(e l; e^ +2 ) mo d n) 
holds. In all four iterations d pre = d post = false holds, resulting in 
two vertices resulting from the intersection of adjacent inequalities: 
{v 1 } = intersect(e i ,e( i+ i) modn ) = intersect(e^ i+ ^ mod „,e( i+2 ) modn), 
{v 2 } = intersect(e( i+2 ) modn,<Z{i+$) modn) = inter sect(e( l+i ) moduli)- 
Hence \E] = conv({v 1 ,v 2 }) + cone({(0, 0)}). 

• Suppose that E is two dimensional and that none of the preceeding cases 
apply. 

— Since c never holds for \E\ > 4 the previous cases deal with all e, e' £ 
E where c(e, e') holds. Because {Ej does not consist of a single point 
described by three inequalities, it must be two dimensional. 

— Let v E V — vert(lEj). By Lemma 13.21 there exists ej £ E such 
that e l Ze (i+1) modn < ir and {v} = intersect^, e (l+1) modn )- Then 
dp OS t is false in iteration i, hence v € P, thus V C P. By Lemma 1X^1 
P Q [£<]• If for all i, e;^e( i+1 ) mod n < n holds then [£] = conv(P) + 
cone((0, 0)). Otherwise let i G [0, m — 1] such that eiZe( i+1 ) mod n > 
it. Then d post will be true in iteration i and d pre will hold in iteration 
(i + 1) mod n yielding rays that are collinear to i\ = (—bi,ai) and 
r-2 = (b(i+i) modn,-a( l+ i) modn)- Since the previous cases do not 
apply, whenever eiZe(i + i) mod n = ""j \E\ > 3 hence Lemma [3.41 can 
be applied. It follows that raydE 1 ]) = cone({ri, T2}), hence [E 1 ] = 
conv(P) + cone{{r\, r 2 })- 



3.3 Reconstruction 

One advantage of the point and ray representation (over one that makes lines 
explicit) is that extreme naturally generates lines as two opposing rays in inde- 
pendent iterations, thus an explicit line case is not required. The remainder of 



the hull function combines the points and rays of the two input polyhcdra to 
construct a corresponding set of inequalities. The advantage of the simplified 
representation carries over to the reconstruction phase in that opposing rays 
from different polyhcdra need not be recognised and reconstituted as a line. 

Theorem 3.2 Given Ei,E 2 C Lin be non-empty, finite, satisfiable and non- 
redundant, E rcs = hull{E\, E 2 ) is non-redundant and the smallest E les C Lin 
with [£ rcs ] D [Si] U IE 2 ] such that for all E G Lin, {Ej D {Eij U \E 2 \ 
lE les }ClE}. 

Proof. The case for Ei — or E 2 = on line 27 is trivial, thus assume 
E 1 U E 2 ^ and that lines 28-35 are executed. Note that for all r G R, \r\ = 1 
due to the normalisation in extreme. Thus (P + {2^/2sr \ r G R})r\(—s, s) 2 = 0. 
Hence after line 12 the predicate inBox(s,p) holds whenever p^P = Q\(P + 
{2y/2sr I r G R}). Line 38 handles the zero dimensional case, hence from line 
39 onwards |Q| > 1. The arithmetic mean q p of all points is calculated in 
line 40. This point serves as reference when comparing two points for counter- 
clockwise ordering. Observe that q p G conv(Q), in particular q p is not on its 
boundary if conv(Q) is two dimensional. The latter ensures for all boundary 
points q\, q 2 G conv(Q), 8(q p , qi) ^ 8(q p , q 2 ), thus line 18 yields a total ordering 
on the boundary points in the two dimensional case. Line 44 performs the classic 
Graham scan which identifies the vertices of conv(Q). 

• To show [£ re J D [Si] U [E 2 j. In particular, show [£ rcs ] D [£1], i.e. for 
allee^ res , [{e}] 2 [Eh]. 

— To show conv(Q) C [{e}]. Suppose e is added in line 60. Then 
the flag add was true and e = connect(qk i , 9fe (i+1) mod m )- For the 
sake of a contradiction, suppose there exists q\~. ^ [{e}] such that 
j £ {i,(i+l) modm}. W.l.o.g. let 0(q p ,q kj ) < 6(q p ,q ki ). There 
exists a line through q p and q kj which intersects the boundary of e at a 
point z. Then q kj G conv({z, qk {i+1) mod m }) which contradicts that q ki 
is a vertex. Since all q ka , ■ ■ ■ , 9fc m _i G [{e}] it follows that conv(Q) C 
[{e}]. Furthermore it is easy to verify that the inequality e added in 
line 58 or 59 holds for both {q^^kA = Q, hence conv(Q) C [{e}]. 

— To show cone(R) C roy([{e}]). 

* Suppose m > 2. If either add was set at line 50 or line 54, there 
exists qj G Q which saturates e such that inBox(s,qj) holds, 
hence q 3 G P. Let r G R. Then qj + 2^sr G conv(Q) C [{e}]. 
Since qj saturates e, it follows that r G ray([{e}]). 

* Suppose m = 2. Assume that inBox(s, q ko ) and inBox(s,q kl ) 
both do not hold, then there exist p G P and r G R with q ko = 
p + 2\[2sr. Note that since m = 2, p saturates e. Continue as 
above. Assume e was added in lines 58-59. Observe that qj = q ki 
saturates e. Again, continue as in the first case. 



Thus [Si] = conu(Pi) + cone(i?i) C conv(P) + cone(R) C [{e}]. Similarly 
for E 2 . Thus lE rcs j D {E^U {E 2 j. 

• To show for all E C Lin, [P] D [Pi] U [P 2 ] [P rc s] C [P]. For the sake 
of a contradiction suppose there exists p G [Pros] such that p £ conv(P) + 
cone(R). Hence for all e G P r0 s, P G [{ e }]- Let G conv(P) + cone(R) 
such that |p — p'\ is minimal. Observe that p' is unique due to convexity. 

— Suppose p' G vert(conv(P) + cone(R)). Hence there exists i G [0,to — 
1] such that p' = qki • Assume that inBox(s, ) does not hold. Then 
there exists p" G P, r G R and A > such that = p" + Ar. 
Since [-Eros] 2 conv(P) + cone(R), p" + 2Ar G [Pros], thus q fei G 
conv({p" ,p" + 2Ar}) which is a contradiction, hence inBox{s,qk i ) 
holds. Thus the flag add is set on line 54 in loop iteration (i — 
1) mod to and i, hence the inequalities e^_ 1} mod „ = qj^~ ) modm ,Qki 
and ei = qk z ,qk {%+1) moi ~ are added to P ros . 

* Suppose ei_iZei < n. Then {q^} = intersect(ei-i,ei). Due to 
convexity |gfe (i _ 1)modm - p| > \qk t ~ P\ and |t?fc (i+1) mwi m - p\ > 
\qki — p\. Hence p' = qk t is the closest point to p in the space 
[{e(i-i) moduli}}. But this implies that p g [P r0 s] which is a 
contradiction. 

* Suppose ej_iZej = it. Then qk (i _ 1} mod m = qk (l+1) modm , thus to = 
2. Note that p G [{e^-i,^}] otherwise p £ P r0 s- Since % i G P, 
inBox(s, q^ ) holds and line 58 or 59 adds an inequality e G P rC s 
with {q^} — inter sect(e,i-i, e) = intersect(e,ei). Observe that 
V £ [{e}] otherwise p G conv({q ki , q fe(i+1) mod m }), thus p £ {E Tes \. 

— Suppose p' £ vert(conv(P) + cone(R)), thus p' is a boundary point of 
conv(P) + cone(R). There exists a line Let 2 such that p' G L and 
conv(P) + cone(R)\L is convex (but not closed). There exists a loop 
iteration i such that the boundary of e = connect(qki , qk {i+1) moim ) is 
exactly L. Since conv(P) + cone(R) C [{e}], p £ [{e}]. 

* Suppose inBox(s, qkj or inBox{s,qk (i+1) mod m ) holds. The flag 
add is true, thus e G P res and p ^ [Pros]- 

* Suppose neither inBox(s,qki) nor inBox(s, qk, i+1)m od m ) holds. 
There exists q,,- G P such that q^ = q rj + Ar for some A > and 
r G -R. Due to the ordering of the points, ki < j < fc(j+i) mo <2 m 
whenever conv(P) + cone(R) is two dimensional. Thus add is set 
in line 54. If conv(P) + cone(R) is one dimensional, to = 2 and 
add is set in line 49. In both cases e G P rC s and thus p £ [{e}]. 

It remains to show that P ros is non-redundant. Note that the loop at lines 
46-62 iterates once for each vertex q^. , creating inequalities = , qk (i+1) moim - 
For all i ^ j, 9(ei) ^ 0(ej) due to the fact that qk m _ 1 are consecutive 

vertices of the hull of conv{Q), thus {eo, . . . , e m _i} has no redundancies. Now 
consider the inequality e\ added at line 58 or 59. Since to = 2, eiZe( i+1 ) moc i m = 



7r and since 6(ei) < 0{e' i ) < 9(e^ + ij m odm), the inequality e\ is not redundant. 
Hence E rcs C {eo, e' , ... , e m _i, e^„_ 1 } is non- redundant. ■ 

The running time of the algorithm is 0(n log n) where n = \Ei | + . After 
the sorting step at line 2 in extreme, each inequality generates at most two rays 
and one point. Thus extreme is in O(nlogn). The flag d pre is true if the angle 
between two consecutive inequalities is at least tt. Thus d pre can only be true 
in at most two loop iterations. Similarly for d post . Hence extreme returns at 
most four rays for each polyhedron, thus \R\ < 8 at line 31 and 0(|Q|) = 0{n). 
The dominating cost in the hull function is the sorting step at line 42 which we 
assume is in 0(n log n). The scan is linear 'A\ and partitions the point set Q into 
vertices and non-vertices. The loop at lines 46-62 runs once for each vertex, 
whereas the inner loop at lines 52-56 runs at most once for each non-vertex. It 
follows that the overall running time is in O(nlogn). 

4 Conclusion 

An 0{n log n) algorithm for calculating the convex hull of planar _ff-polyhedra 
has been presented, thereby improving on existing approaches. The algorithm 
applies a novel box construction which reduces the problem to calculating the 
convex hull of a set of points. Implementing the algorithm exposed a number 
of subtleties which motivated a complete proof of the algorithm. 
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